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ABSTRACT 

XMM-Newton X-ray spectra of the hard state Black Hole X-Ray Binaries (BHXRBs) 
SWIFT J1753. 5-0127 and GX 339-4 show evidence for accretion disc blackbody emis- 
sion, in addition to hard power-laws. The soft and hard band Power-Spectral Densities 
(PSDs) of these sources demonstrate variability over a wide range of time-scales. How- 
ever, on time-scales of tens of seconds, corresponding to the putative low-frequency 
Lorentzian in the PSD, there is additional power in the soft band. To interpret this 
behaviour, we introduce a new spectral analysis technique, the 'covariance spectrum', 
to disentangle the contribution of the X-ray spectral components to variations on dif- 
ferent time-scales. We use this technique to show that the disc blackbody component 
varies on all time-scales, but varies more, relative to the power-law, on longer time- 
scales. This behaviour explains the additional long-term variability seen in the soft 
band. Comparison of the blackbody and iron line normalisations seen in the covari- 
ance spectra in GX 339-4 implies that the short-term blackbody variations are driven 
by thermal reprocessing of the power-law continuum absorbed by the disc. However, 
since the amplitude of variable reflection is the same on long and short time-scales, we 
rule out reprocessing as the cause of the enhanced disc variability on long time-scales. 
Therefore we conclude that the long-time-scale blackbody variations are caused by in- 
stabilities in the disc itself, in contrast to the stable discs seen in BHXRB soft states. 
Our results provide the first observational evidence that the low-frequency Lorentzian 
feature present in the PSD is produced by the accretion disc. 
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1 INTRODUCTION 

X-ray spectra of BHXRBs show evidence for a two-phase 
structure to the accretion flow, an optically thick, geo met- 
rically thin accretion disc (|Shakura fc Sunvaevl I1973T ) giv- 
ing rise to a blackbody component in the X-ray spectrum, 
and a hot optically-thin component, modelled as a power 
law. The relative strengths of these two co mponents de- 
fine t he appearance of different 'states' (e.g. iHoman et al.l 
l200ll h In the hard state, which we focus on in this paper, 
the power-law emission dominates the total luminosity. It 
has been suggested that the power-law is produced by an 
inner, opticall y-thin Advection Dominated Accreti on Flow 
(ADAF) (e.g. lEsin. McClintock. fc Naravanl 1 19971 ) . which 
replaces the inner optically thick disc at low accretion rates, 
and extends down to the Innermost Stable Circular Or- 
bit (ISCO). This implies that the optically thick disc is 
truncated at some transition radius. Alternatively, the op- 
tically thick and optically thin components may co-exist 
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over some range of radii, e.g . if the thin disc is sand- 
wiched by a hot flow o r corona jWitt. Czernv fc Zvcki|[l997l ; 
IChurazov et al.l 1200 il l . The corona may in turn evaporate 
the innermost regions of the disc at low accretion rates ( e.g. 
iMever-Hofmeister fc Meverl2003l ; lMaver fc Pringlel2007h . so 
that the optically thin flow is radially separated from the op- 
tically thic k disc as for the ADAF model. It has also been 
sugg ested (|Malzadl2007h that a cold thin accretion disc ex- 
tends close to the black hole in the hard-state. In this model, 
most of the accretion power is transported away from the 
disc to power a strong outflowing corona and jet. Since the 
ADAF and corona perform similar roles in that they upscat- 
ter soft photons to produce the observed power-law, we shall 
henceforth refer to both components interchangeably as the 
corona, without necessarily favouring either picture. 

In the hard state, the variability of the power- 
law continuum offers further clues to the structure of 
the accretion flow. Studies of the timing properties of 
hard state BHXRBs show that the frequencies of broad 
Lorentzian features in their Power-Spectral Density func- 
tions (PSDs) correlate with the strength of the reflec- 
tion features as well as the steepness of the power 
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law continuum llRevnivtsev. Gilfanov"fc~ Churazovl l200ll; 
iGilfanov. Churazov fc Revnivtsevll 19991 ). These correlations 
can naturally be explained if the Lorentzian frequencies cor- 
respond to a characteristic timescale at the disc truncation 
radius, e.g. the viscous time-scale, so that as the truncation 
radius increases the Lorentzian frequency decreases, with 
disc reflection and Compton cooling of the optically thin 
hot flow by disc photons decreasing accordingly. In this pic- 
ture, the truncation radius of the thin disc acts to generate 
the signals of the lowest- frequency Lorentzian in the PSD, 
while the highest-frequency Lorentzians may be generated 
at the inner most radius of the hot inner coronal flow, i.e., 
at the ISCO (|Done et alj|2007h . 

Regardless of whether the corona is radially or verti- 
cally separated from the thin disc, photons upscattered by 
the corona should interact with the disc. This interaction 
gives rise to reflection features in the X-ray spectrum, in- 
cluding fluorescent iron line emission and a reflection con- 
tinuum due to Compton scattering off the disc material. An 
often-neglected consideration is that a significant fraction 
of the photons interacting with the disc are absorbed and 
the disc heated in a process known as thermal reprocessing. 
Provided that the disc subtends at least a moderate solid 
angle as seen by the corona, this effect should be particu- 
larly significant in the hard state, where the coronal power 
law continuum dominates the total luminosity. When the 
power-law luminosity impinging on the disc is high com- 
pared to the disc luminosity due to internal heating, then a 
significant fraction of the disc blackbody emission should be 
reprocessed and will therefore track variations of the power 
law continuum. 

The anticipated correlated variations of the blackbody 
and power-law emission can be studied using variability 
spectra, e.g. the rms spectrum, w hich show only the v ari- 
able components of the spectrum (jVaughan et al. 1 120031 ). If 
the geometry is such that the observed power-law produces 
reprocessed blackbody emission by X-ray heating the disc, 
then both power-law and blackbody components should ap- 
pear together in the variability spectra. Furthermore, by 
selecting different time ranges covered by these variability 
spectra (i.e., analogous to the method of Fourier-re solved 
spectrosopy. iRevnivtsev. Gilfanov fc Churazovl Il999), it is 
possible to determine whether the low-frequency part of 
the PSD has a different origin to the high-frequency part 
in terms of the contributions of blackbody and power-law 
components. If the optically thick disc does drive the low- 
frequency Lorentzian, we predict that the blackbody com- 
ponent should be stronger in the corresponding variability 
spectrum. 

In this work, we examine the variability spectra of two 
hard state black hole X-ray binaries, SWIFT J1753. 5-0127 
and GX 339-4, which have good XMM-Newton data for 
studying simultaneous variations of disc and power-law. Pre- 
vious analyses of time-averaged spectra for these XMM- 
Newton observations have shown the existence of the ex- 
pected blackbody components, together with relativistically 
broadened reflection features, which have been used to ar- 
gue that the disc is truncated at significantly smaller radii 
than previously thought , perhaps only a few gravitationa l 
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at larger radii. In the following Section, we describe the ob- 
servations and data reduction. In Section[3]we show the soft 
and hard-band PSDs obtained from the data, and present 
a technique to produce a type of rms spectrum, the 'covari- 
ance spectrum' which we use to identify the variable spec- 
tral components for each hard state source. In particular 
we show that, although both power- law and disc blackbody 
emission are correlated, as expected from thermal reprocess- 
ing, the disc is relatively more variable than the power-law 
on longer time-scales (corresponding to the low-frequency 
Lorentzian), contrary to what we would expect from simple 
reprocessing models. In Section [3] we discuss the interpreta- 
tion of our results and present further analysis to suggest 
that the low-frequency Lorentzian corresponds to fluctua- 
tions intrinsic to the disc. A summary of our conclusions is 
given in Section [5] 



2 OBSERVATIONS AND DATA REDUCTION 

SWIFT J1753. 5-0127 was observed during revolution 1152 
for 42 ks by XMM-Newton EPIC-pn on 2006 March 24 
in pn-timing mode using the medium optical filter. The 
events list was screened using the perl script xmmcleaivQ 
to select only events with FLAG=0 and PATTERN sC 4. 
Examination of light curves showed no evidence for back- 
ground flaring. Using SAS Version 7.0, evselect was used 
to extract the mean spectrum, f ollowing the procedure of 
iMiller. Homan. fc Miniuttil (|2006l ). using event positions be- 
tween 20 and 56 in RAWX and using the full RAWY range. 
The SAS tasks arfgen and rmfgen were used to gener- 
ate the Ancillary Response File (ARF) and Redistribution 
Matrix File (RMF). 

GX 339-4 was observed by XMM-Newton on 2004 
March 16 during revolutions 782 and 783, again in pn- 
timing mode using the medium optical filter. The SAS com- 
mand evselect was used to filter the events on TIME to 
avoid background flaring, producing a total combined expo- 
sure of 127 ks. Data reduction was very similar to that for 
SWIFT J1753.5-0127, using SAS version 7.0, but care was 
taken to avoid pile-up. Successive columns from the centre 
of the image were excised in RAWX and spectra extracted 
until no discernable difference in spectral shape could be 
identified between successive selections on RAWX, implying 
that pile-up is no longer significant. The EPIC-pn data was 
deemed free of pile-up using extraction regions in RAWX 
from columns 30 to 36 and 40 to 46 (i.e., columns 37 to 39 
inclusive were excised). Background spectra were selected in 
RAWX from columns 10 to 18 over the full RAWY range. 
To generate an appropriate ARF for the data made in this 
way, it was necessary to use arfgen to generate an ARF for 
the full region in RAWX from columns 30 to 46 and gener- 
ate a second ARF from the spectrum of the excluded region 
(RAWX columns 37 to 39) and then subtract the latter from 
the former using the command ADDARlJf). 

For both sources, the FTOOL GRPPHA was used to scale 
the background spectra and to ensure that a minimum of 20 
counts were in each bin for x 2 fitting. The EPIC-pn covers 
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the energy range from 0.2-10 keV, but the EPIC Calibration 
Status Document^ recommends restricting the fit to energies 
greater than 0.5 keV. To be conservative, our EPIC-pn spec- 
tral fits were restricted to the range 0.7 to 10.0 keV. 

We extracted data from RXTE observations which were 
contemporaneous with the XMM-Newton observations. We 
used the high time and spectral-resolution PCA Event mode 
data to extract mean spectra and make rms and covariance 
spectra using the same time binning as the EPIC-pn data 
(see Section EH). Throughout this work we use RXTE PCA 
data in the 3-25 keV range, and fit these together with the 
corresponding EPIC-pn spectra, tying all fit parameters to- 
gether but allowing RXTE spectral fits to be offset by a 
constant factor with respect to the simultaneous EPIC-pn 
fits, to allow for the difference in flux calibration between 
the PCA and EPIC-pn instruments, and also for the fact 
that slight flux differences may result from the fact that 
the RXTE observations covered shorter intervals than the 
XMM-Newton data. A 1 per cent systematic error was as- 
sumed in all spectral fits to account for uncertainties in 
instrumental response and cross-calibration. The signal to 
noise of RXTE HEXTE data was not sufficient to perform 
the covariance analysis we describe in section 3.2. 



3 ANALYSIS AND RESULTS 
3.1 Power-spectral densities 

Before describing the spectral analysis method and results, 
we first examine the timing properties of each source in soft 
and hard spectral bands using the power-spectral density 
function (PSD). For each source, we used the EPIC-pn data 
to generate light curves with 1 ms time binning in two en- 
ergy bands: 0.5-1 keV (soft) and 2-10 keV (hard). Only com- 
plete segments of 131 s duratiorQ were used to construct the 
PSDs, segments with gaps (which can be common at high 
count rates in timing mode) were skipped over. The result- 
ing Poisson-noise-subtracted PSDs are shown in figure [1] It 
is clear from the figure that the PSDs are significantly differ- 
ent in shape between the two bands. In GX 339-4, the hard- 
band PSD shows two broad components, with the higher- 
frequency component appearing to be shifted to even higher 
frequencies in the hard band. In SWIFT J1753. 5-0127, soft 
and hard PSDs appear, within the noise, to overlap more 
closely at higher frequencies. However, in both sources the 
soft band PSD shows relatively larger low-frequency power 
compared to the PSD components at frequencies above 
~ 0.2 Hz. This ~ 20% extra low-frequency power at lower 
X-ray energies may be associated with the soft excess emis- 
sion, e.g. if the disc is varying more than the power-law on 
longer time-scales. Alternatively, there could be power-law 
spectral variability that applies only on longer time-scales 
and causes the soft band variability amplitude to be en- 
hanced, i.e., due to steepening of the power-law spectral 
slope on longer time-scales. To determine the origin of the 
extra low-frequency power in the soft band, we must carry 



3 http: / /xmm2.esac.esa.int/external/xmm_sw_cal/calib/index.shtml 

4 Specifically, containing 131072 time bins. 



out a spectral analysis of the variations on different time- 
scales, the methodology for which we describe in the next 
section. 

3.2 Methodology: rms and covariance spectra 

When fitting models to X-ray spectra, it is typical to use 
mean X-ray spectra that only describe the time-averaged 
spectral shape of a source. Such fits say nothing about the 
way the different spectral components such as disc black- 
body and power-law vary with respect to each other in time. 
By looking at the absolute amplitude of variations in count 
rate as a function of energy we can construct 'variability 
spectra' which pick out only the time-varying components. 
One technique is to construct Fourier- frequency resolved 
spectra, obtaining a PSD for each individual energy chan- 
nel and integrating the PSD over a given Fourier-frequency 
range in order to measure the variance in that channel, which 
is used to obtain the rms and so construct the spectrum due 
to components which vary over that fr equency range (e.g. 
iRevnivtsev. Gilfanov fc Churazovl I1999T ). This approach is 
attractive in that it allows the user to pick out complex pat- 
terns of spectral variability where components have different 
time-scales of variation, as may be implied by the soft and 
hard band PSDs shown in figureQ] Here we will consider only 
two time-scales of variability, corresponding roughly to the 
high and low-frequency parts of the PSD which show signifi- 
cant relative differences between the soft and hard bands. To 
approximate the more complex Fourier-resolved approach 
we will measure over two time-scale ranges a variant of the 
'rms spectrum', which measures the absolute Root-Mean- 
Squared variability (rms) as a function of energy (e.g. see 
IVaughan et"afll2003h . 

Producing the rms spectrum involves allocating each 
photon event to a time and energy bin, dividing the light 
curve into segments consisting of TV time bins per segment 
and then working out the variance in each segment for each 
energy bin according to the following standard formula: 

JV 

where Xi is the count rate in the i th bin and X is the mean 
count rate in the segment. The expectation of Poisson noise 
variance, given by the average squared-error crf lr) is sub- 
tracted, leaving the 'excess variance', a^ s , in each segment. 
The excess variances can then be averaged over all of the 
segments of the light curve. The square root of the aver- 
age excess variance plotted against energy forms the rms 
spectrurrQ. By selecting the time bin size and the segment 
size, we can isolate different time-scales of variability, effec- 
tively replicating the Fourier resolved approach for the two 
time-scale ranges that we are interested in. For this purpose, 
we choose two combinations of bin size and segment size. To 
look at variations on shorter time-scales we choose 0.1 s time 
bins measured in segments of 4 s (i.e., 40 bins long), i.e., cov- 
ering the frequency range 0.25-5 Hqj. For longer time-scale 

5 This is absolute rms, since we are not normalising by the mean 
count rate. 

6 The upper frequency limit is set by the Nyquist frequency, 
^Nyq = 1/(2 At) where At is the time bin size. 
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Figure 1. Soft and hard band PSDs for SWIFT J1753. 5-0127 and GX 339-4. The solid line denotes the 0.5-1 keV PSD, while the dotted 
line denotes the 2-10 keV PSD. The Poisson noise component has been subtracted from all PSDs. 



regions we use 2.7 s bins in segments of 270 s, i.e., covering 
the range 0.0037-0.185 Hz. Note that these two frequency 
ranges do not overlap and also cover the two parts of both 
source PSDs which show distinct behaviour in soft and hard 
bands. 

There is, however, a problem with the rms spectrum 
when signal to noise is low, e.g. at higher energies. It is 
possible for the expectation value of the Poisson variance 
term to be larger than the measured average variance term, 
producing negative average excess variances. If this is the 
case, it is not possible to calculate the rms at these energies 
and this introduces a bias towards the statistically higher- 
than-average realisations of rms values, which can still be 
recorded. In order to overcome these problems we have de- 
veloped a technique called the 'covariance spectrum'. The 
covariance is calculated according to the formula: 



N 



— J2(Xi-X)(Yi-Y) 



(2) 



where Yi now refers to the light curve for a 'reference band' 
running over some energy range where the variability signal- 
to-noise is large. In this work, we use reference bands of 1-4 
keV for EPIC-pn data and 3-5 keV for PCA data. In other 
words, the covariance spectrum is to the rms spectrum what 
the cross-correlation function of a time series is to its auto- 
correlation function. The covariance spectrum therefore does 
not suffer from the same problems as the rms spectrum, as 
no Poisson error term has to be subtracted, since uncorre- 
cted noise tends to cancel out and any negative residuals 
do not affect the calculation. To remove the reference band 
component of the covariance, and produce a spectrum in 
count-rate units, we obtain the normalised covariance for 
each channel using: 



' cov.norm 



2 

xs,y 



(3) 



where <j xs „ is the excess variance of the reference band. 
Therefore, the only requirement for there being a valid, un- 
biased value of covariance at a given energy is that the refer- 
ence excess variance is not negative. This is usually the case, 
since the reference band is chosen to include those energies 



with the largest absolute variability. When the covariance is 
being calculated for an energy channel inside the reference 
band, the channel of interest is removed from the reference 
band. The reasoning behind this is that if the channel of in- 
terest is duplicated in the reference band, the Poisson error 
contribution for that channel will not cancel and will con- 
taminate the covariance. One can think of the covariance 
technique as applying a matched filter to the data, where the 
variations in the good signal-to-noise reference band pick out 
much weaker correlated variations in the energy channel of 
interest that are buried in noise. In this way the covariance 
spectrum picks out the components of the energy channel of 
interest that are correlated with those in the reference band. 

It is important to note that the covariance spectrum 
only picks out the correlated variability component and is 
therefore a more appropriate measure than the rms spec- 
trum in constraining the reprocessing of hard photons to 
soft photons, which will result in correlated variations. When 
the raw counts rms and covariance spectra are overlaid, as 
in figure [5] they match closely indicating that the reference 
band is well correlated with all other energies (the sp ectral 
'coherence' is high, e.g. see IVaughan fc Nowakl U~997l . An- 
other advantage of the matched filter aspect of the covari- 
ance spectrum is that it leads to smaller statistical errors 
than the rms spectrum. Specifically, the errors are given by: 



Err [a 



'xs,x u crr.y 



■,x + a\ 



crr.x '-'crr.y 



cov .norm 



NMai a , y 



where M denotes the number of segments and subscripts x 
and y identify excess variances and Poisson variance terms 
for the channel of interest and reference band respectively. 
This error equation can be derived simply from the Bartlett 
formula for the error on the zero-lag cross-correlation func- 
tion, assuming th at the source li ght curves have unit y in- 
trinsic coherence l|Bartlett|[l955l: iBox fc Jenkins! 1 197rJ) . By 
comparison with equation B2 of IVaughan et all 20031 . using 
the relation: 



<< 



(which is true because the reference band has good signal to 
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Figure 2. Overlaid covariance and rms spectra for EPIC-pn and 
PCA data for SWIFT J1753.5-0127. The large error bars in the 
EPIC-pn data indicate negative excess variances where no rms 
value could be determined. 

noise) it can be shown that the errors on the covariance are 
smaller than corresponding errors on the rms values. 

Finally, we note that, since the covariance is analogous 
to the zero-lag cross-correlation function, it could be affected 
by intrinsic time lags in the data. Using measurements of the 
cross-spectral phase-lags between various energy bands, we 
have confirmed that the lags between hard and soft band 
variations are smaller than the time bin sizes used to make 
the long and short time-scale covariance spectra. Thus, in- 
trinsic time lags will have no effect on our results. 

3.3 SWIFT J1753. 5-0127 

We first consider the covar iance spectra of SWIFT J 1753. 5- 
0127, and use xspec vlsQ (|Arnaudlll996t ) to fit the long and 
short-time-scale data together with the mean spectrum in 
order to identify the origin of the additional long-time-scale 
variability which can be seen in the soft band PSD. The co- 
variance spectra do not show sufficient signal-to-noise to de- 
tect the rathe r weak iron line present in the mean sp ectrum 
of this source (jHiemstra et al.ll2009l ; iReis et alj |2009) , so for 
simplicity we fit only a simple power-law and multicolour 
disc blackbody diskbb, together with neutral absorption. 
We fit the short and long-time-scale covariance and mean 
spectra simultaneously, tying the absorbing column density 
and renormalising constant for PCA data to be the same for 
the mean and covariance spectra. An F-test showed that the 
disk blackbody temperature does not change significantly 
between spectra, so that was also tied to be the same for 
all spectra. The remaining parameters were allowed to be 
free between the covariance and mean spectra. For those 
parameters that were free to vary, EPIC-pn and PCA val- 
ues were tied together. The \ 2 °f the final fit was 1802 for 
1948 degrees of freedom (d.o.f.) and full fit parameters are 
listed in Table [1] The best-fitting unfolded spectra and the 
corresponding data/model ratios are shown in figure [3] 

7 http:/ /heasarc. gsfc.nasa.gov/xanadu/xspec/ 



Table 1. Fit parameters for mean and covariance 
spectra of SWIFT J1753. 5-0127 for the model CON- 
stant*phabs*(diskbb+powerlaw) 



Parameter 


Mean 


2.7s-270s Cov 


0.1s-4s Cov 


Cpca 

r 

A:Tj n 

Adisc 


-, r, 99 +o.a06 
r.zzz_ 00g 

194+0 007 

1 «oc:+0-005 

n n7"*i +0-0009 
u.u(oi_ 0009 

0-286±S;g« 

320 ± 25 


tied to mean 
tied to mean 
1.55 ±0.02 

U-UJ-J.8_o.oo04 
tied to mean 

129 ±11 


tied to mean 
tied to mean 
1.55 ±0.03 

n ni98+° ooo 6 
U.Ulz8_ 00Q6 

tied to mean 

75 ±20 



From top to bottom, the parameters are: constant renormalising 
factor applied to model fit to PCA data; neutral absorbing 
column density (10 22 cm -2 ); power-law photon index; 
power-law normalisation (photons cm -2 s _1 keV~ 1 at 1 keV); 
disc blackbody inner radius temperature (keV); disc blackbody 
normalisation ([r; n /(Z)/10 kpc)] 2 where n n is the disc inner 
radius in km and D is the distance in kpc). All errors are 90 per 
cent confidence limits. Since disc temperature and normalisation 
are highly correlated, the errors on blackbody normalisation are 
obtained while fixing temperature at the best-fitting value. 

To interpret the fits to the covariance spectra, consider 
the case where the observed PSDs have identical shapes in 
both hard and soft bands (they may have different normal- 
isations). In this case, the ratio of soft to hard-band vari- 
ability amplitudes measured over the same time-scale range 
will be identical for any given time-scale range. Therefore, 
since (for high-coherence variations) the covariance spec- 
trum quantifies variability amplitude as a function of energy, 
the shape of the covariance spectrum will be independent 
of time-scale. On the other hand, if the soft band contains 
more long-time-scale variability relative to short-time-scale 
variability than the hard band, the long-time-scale covari- 
ance spectrum will appear softer than the short-time-scale 
covariance spectrum: lower energies show correspondingly 
greater variability on long time-scales and therefore larger 
fluxes in the covariance spectra. Just such an effect is seen 
in the model fits to the covariance spectra: the long-time- 
scale covariance spectrum is softer than the short-time-scale 
covariance spectrum, because the disc blackbody normali- 
sation is higher on long time-scales. The power-law slope is 
remarkably similar on both long and short time-scales how- 
ever. Therefore the additional long-term variability in the 
soft-band PSD seems to result from additional variability of 
the disc blackbody, not any extra power-law variability (e.g. 
due to spectral pivoting). 

We can show in a model-independent way that the 
blackbody component is the source of the additional long- 
term variability seen in the soft-band PSD, by measuring the 
'covariance ratio', i.e., plotting the ratio of long-time-scale 
to short time-scale covariance spectra. By plotting the ratio, 
we remove most of the effects of the spectral response (ex- 
cept the smearing by the instrumental resolution) and can 
see more clearly the nature of any spectral differences. The 
covariance ratio plot is shown in figure [4] and shows clearly 
that the excess in long-time-scale covariance arises abruptly 
below 2 keV, as expected if caused only by the disc emis- 
sion, rather than being due to a gradual change across the 
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Figure 3. Top panel: Unfolded spectra for SWIFT J1753.5- 
0127 showing the relative contributions of the disc blackbody 
and power law components in the mean spectrum and short and 
long-time-scale covariance spectra. Dotted lines indicate spectral 
components of the short time-scale covariance, dashed lines in- 
dicate those of the long time-scale covariance. Bottom panel: 
SWIFT J1753. 5-0127 data/model ratios for the best-fitting disc 
blackbody plus power-law model. 



whole spectrum which might be caused by differences in the 
power-law index. The interpretation of the short time-scale 
covariance spectrum is addressed later in section 4.4. 

Having established that the reference band is well cor- 
related with all other energies, the normalisations of the co- 
variance spectral components, given in Table Q] relative to 
those of the mean spectrum indicate the rms variability am- 
plitude. For example, dividing the disc blackbody covariance 
normalisations by that of the mean spectrum indicates frac- 
tional rms for the disc emission of 40 per cent and 23 per 
cent over long and short time-scale ranges respectively. The 
power-law variability is more complicated to quantify, since 
the covariance spectra show significantly harder power-law 
slopes than the mean. This result implies that the power-law 
fractional rms increases with energy. The power-law normal- 
isations imply that at 1 keV the fractional rms of the power- 




• XMM Cov Ratio 2. 7s/ 0.1s 
□ XTE Cov Ratio 2.7s/0.1s 
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Figure 4. Covariance ratios for SWIFT J1753. 5-0127 show- 
ing the extra blackbody component which contributes on longer 
timescales. 



law is 16 per cent and 18 per cent over the long and short 
time-scale ranges respectively. 



3.4 GX 339-4 

We next produced mean, long and short time-scale covari- 
ance spectra for EPIC-pn timing mode and PCA observa- 
tions of GX 339-4. We carried out fits to the spectra using a 
model consisting of a narrow Gaussian emission line (width 
a = 0.01 keV) and power law (including the multiplicative 
model hrefl to account for the reflected continuum from 
a cold optically thick disk), smeared with the convolution 
model kdblur in order to account for relativistic broaden- 
ing of the reflected emission fr om material orbi ting a maxi- 
mally spinning black hole (see lReis et al.ll2008l for evidence 
for a high black hole spin in GX 339-4). A multicolour disc 
model diskbb was also included to account for the soft X- 
ray excess. Strictly speaking this component should also be 
smeared, but it is not appropriate to smear this model with 
kdblur, as the disc temperature is itself a function of ra- 
dius. Since in the kdblur model the disc innermost radius 
and iron line radial emissivity are somewhat degenerate, the 
emissivity index was fixed at -3 (i.e., disc-like). Also, the 
outer radius was fixed at 400 Rg- The remaining param- 
eters, except the black body normalisation, reflection cov- 
ering fraction, photon index, power law normalisation and 
Gaussian normalisation were tied to be the same between 
the mean and both covariance spectra (we confirmed with 
F-tests that the tied parameters do not vary significantly 
between the spectra). As with SWIFT J1753. 5-0127, the 
EPIC-pn and PCA parameters were tied together for the 
same spectra and a renormalising constant factor was in- 
cluded to allow for differences in EPIC-pn and PCA flux 
calibration. The final fit had a X 2 of 2489 for 2199 d.o.f. and 
a full list of model parameters is given in Table [2] The best- 
fitting unfolded spectra and data/model ratios can be seen 
in figure [S] 

The fits to the covariance spectra for GX 339-4 show 
a similar pattern to that seen for SWIFT J1753.5-0127, in 
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Figure 5. Top panel: Unfolded spectra fitted for GX 339-4 show- 
ing the relative contributions of the blackbody, power law (incor- 
porating reflection continuum) and iron line components. Dotted 
lines indicate spectral components of the short time-scale covari- 
ance, dashed lines indicate those of the long time-scale covariance. 
Bottom panel: GX 339-4 data/model ratios for the best-fitting 
disc blackbody plus power-law and relativistically smeared reflec- 
tion model. 



that the long time-scale covariance spectrum is softer than 
the short time-scale covariance spectrum because of a sig- 
nificantly stronger disc blackbody component, while their 
power-law indices are very similar. We can confirm this in- 
terpretation using a plot of the covariance ratio, which is 
shown in figure [6] and shows a similar rise at low ener- 
gies to that seen for SWIFT J1753. 5-0127, which underlines 
the interpretation that the disc blackbody component is the 
main contributor to the additional variability seen at low 
frequencies in the soft band in GX 339-4. The reader is re- 
ferred to sections 4.1 and 4.4 for a consistent interpretation 
of the short time-scale covariance spectrum in GX 339-4 
and SWIFT J1753.5-0127. Note that the apparent emission 
feature around 2 keV is probably an instrumental effect, 
possibly related to mild pile-up in timing mode, since it dis- 
appears when larger regions of RAWX are excised from the 



Table 2. Fit parameters for mean and covariance spectra of GX 
339-4 for the model constant * phabs * (diskbb+kdblur(hrefl 
* powerlaw+gaussian)) 



Parameter 


Mean 


2. / S-2 / US bov 


U.1S-4S Uov 


Cpca 


1.296 ±0.003 


tied to mean 


tied to mean 




0.574 ± 0.006 


tied to mean 


tied to mean 


r 


1.766 ± 0.006 


1.896 ±0.010 


1.875 ±0.011 




0.266 ± 0.002 


0.110 ±0.001 


0.107 ±0.001 


kT in 


0.177 ±0.001 


tied to mean 


tied to mean 


A disc /10 4 


8.20 ±0.65 


3.20 ±0.26 


2.46 ±0.21 


Rin 


4.08 ± 0.16 


tied to mean 


tied to mean 




40.2 ± 1.4 


tied to mean 


tied to mean 


£ Fo 


6.42 ± 0.07 


tied to mean 


tied to mean 


A Pc /10- 3 


1.17 ±0.07 


0.61 ±0.07 


0.58 ±0.09 


CF 


0.91 ±0.04 


1.13 ±0.08 


1.16 ±0.09 



For the definitions of the parameters listed in the first six rows 
see Table [T] The additional parameters shown are for the 
relativistically smeared reflection (from top to bottom): disc 
innermost radius (units of Rq)\ disc inclination (fixed to be the 
same in both KDBLUR and hrefl); Gaussian line energy; 
Gaussian normalisation (photons cm~ 2 s" 1 ); covering fraction 
of the reflection (where 1.0 corresponds to 2-7T steradians). 



data. Due to the 1 per cent systematic included in the spec- 
tral fitting, this feature has no effect on the fit results. 

Based on the ratio of component normalisations to those 
in the mean spectrum, the fractional rms values for the disc 
component are 39 per cent and 30 per cent over the long 
and short time-scale ranges respectively. The corresponding 
power-law fractional rms values at 1 keV are 41 per cent and 
40 per cent. Using the same approach, we can also define a 
fractional rms for the iron line emission for long and short 
time-scales, at ~ 50 per cent in each case. Interestingly, the 
~ 25 per cent increase in iron line rms over that of the power- 
law emission (which drives the line variability) is comparable 
to the increase in reflection covering fraction from the mean 
to covariance spectra. These results may imply the pres- 
ence of an additional constant power-law component which 
dilutes the power-law fractional variability but does not con- 
tribute to reflection. It is also interesting to note that the 
power-law component in the covariance spectra for GX 339- 
4 is softer than in the mean spectrum, i.e., the difference is 
in the opposite sense to that seen in SWIFT J1753. 5-0127. 



4 DISCUSSION 

We have shown that disc blackbody emission contributes sig- 
nificantly to the X-ray variability spectra in the hard state of 
the black hole candidates SWIFT J1753.5-0127 and GX 339- 
4, and moreover, that the disc emission is the origin of the 
additional soft band variability seen on longer time-scales 
in both sources, which manifests itself as an enhanced low- 
frequency component in the PSD. In this section, we discuss 
the evidence for a connection between disc and power-law 
variability through X-ray reprocessing, and then consider 
two possible explanations for the enhanced disc variability 
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Figure 6. Covariance ratios for GX 339-4 showing the extra ther- 
mal component on longer timescales 



on long time-scales, in terms of geometry changes or fluctua- 
tions intrinsic to the accretion disc. Finally we will compare 
our results and interpretation to the wider picture of differ- 
ent accretion states. 



4.1 Power-law variability and evidence for 
thermal reprocessing 

The power-law emission clearly dominates the X-ray lumi- 
nosity in both sources, as can be seen in figures [3] and [5] In 
this situation, if the disc sees a reasonable fraction of the 
power-law emission, we should expect that X-ray heating of 
the disc, i.e., thermal reprocessing of power-law emission, 
will produce a significant fraction of the observed disc lumi- 
nosity. In fact, the thermal reprocessed emission is directly 
related to the disc reflection component in the spectrum: if 
a fraction of incident power-law luminosity / is reflected by 
the disc (both through Compton reflection and emission line 
fluorescence), then a fraction 1 — / must be absorbed and 
will be reprocessed into thermal blackbody radiation. The 
reflected fraction / depends on disc ionisation state but sim- 
ple exploration of the disc reflection models in xspec shows 
that it is typically 30-40 per cent of the incident luminos- 
ity, so that around 60-70 per cent of the incident luminos- 
ity is reprocessed into disc blackbody emission. For typical 
iron Ka line equivalent widths (around 1 keV with respect 
to the reflection continuum), we then expect the line flux 
to be of order 1 per cent of the thermally reprocessed flux. 
We might reasonably assume that the blackbody component 
in the short time-scale covariance spectrum is produced by 
thermal reprocessing of the varying power-law, which also 
drives the line emission in the same spectrum. The unab- 
sorbed disc flux is ~ 4.5 x 10 -10 erg cm -2 s _1 , and the iron 
line flux is ~ 6 x 10 -12 erg cm -2 s _1 which is 1.3 per cent 
of the reprocessed flux, i.e., consistent with a reprocessing 
origin for the thermal emission, at least on short time-scales. 

In SWIFT J1753.5-0127, the disc blackbody emission 
is considerably weaker than in GX 339-4, and correspond- 
ingly we would expect a relatively weak iron line, with a 
few tens of eV equivalent width, which is only just con- 



sistent with th e lower -limi ts on line s t rengt hs reported by 
Iffiemstra et all (|2009h and iReis etafl <|2009h . It is possible 
that the disc in SWIFT J 1753.5-0127 is substantially ionised 
(e.g. see lReis et all l2009) which would enhance line emission 
relative to the absorbed (and hence thermally reprocessed) 
emission. 

Since it is likely that there is substantial X-ray heat- 
ing of the disc, one must interpret the observed blackbody 
normalisations in the mean spectra with cautiorjj. The disc 
emissivity may be more centrally concentrated than the the- 
oretically expected R~ 3 law, and so the normalisations indi- 
cate better the emitting surface area and cannot be simply 
translated to an inner radius. Nonetheless, the inferred emit- 
ting areas are still relatively small, implying distance scales 
of tens to hundreds of km (assuming distances > 7 kpc and 
6-15 k pc for SWIFT J1753. 5-0127 and GX 339-4 respec- 
tively; [Zuritaelal]|200i; |HxneseFaL| [iooi) . We also note 
here that although the need for disc blackbody emission to 
explain t he spectrum of SWIFT J1753. 5-0127 has been ques- 
tioned bv lHiemstra et all (|2009l ). the model-independent co- 
variance ratio plots in figure [4] show that a distinct soft com- 
ponent must be present in order to explain the difference in 
the shapes of the covariance spectra. 

It is interesting to note that the power-law indices of the 
covariance spectra are different to those of the mean spectra, 
but in an opposite sense for each of the two sources consid- 
ered here: compared to the mean spectrum SWIFT J1753.5- 
0127 shows a harder power-law in the covariance spec- 
tra, while GX 339-4 shows a softer power-law. The differ- 
ence may be caused by flux-dependent spectral pivoting 
or steepening, so that as flux increases the spectrum gets 
softer in GX 339-4, increasing the covariance at soft ener- 
gies relative to the mean, while the opposite effect occurs in 
SWIFT J1753. 5-0127 (it hardens as it gets brighter). The 
difference in behaviour may be related to the source lumi- 
nosity: if they lie at similar distances SWIFT J1753. 5-0127 
is at least a factor 10 less luminous than GX 339-4, imply- 
ing a significantly lower accretion rate. Correlations between 
flux and spectral-hardness have been seen in BHXRB hard 
states, and interestingly the sign of the correlation appears 
to switch over from negative to positive at low lu minosi- 
ties, both on short time-scales llAxelsson et al.l 20081 ) and in 
the long-term global correlation ( Wu fc Gull2008l ). The same 
switch in flux-hardness correlation could be related to the 
different power-law behaviour of SWIFT J1753. 5-0127 and 
GX 339 which we see here. 

The same pattern appears on both long and short time- 
scales, which show almost identical power-law indices in 
their covariance spectra, so that the effect of the power- 
law spectral variability will be to change the normalisa- 
tion of the PSD, but not the shape. Since most of the 
power-law luminosity will be found at tens of keV (assum- 
ing a thermal cutoff at around 100 keV), spectral steepen- 
ing with flux will cause the observed GX 339-4 0.5-10 keV 
variability amplitude to be enhanced compared to the true 
luminosity variations. Conversely, observed variations in 
SWIFT J1753. 5-0127 will be smaller than the total luminos- 



8 The absolute normalisation in the covariance spectra is a func- 
tion of variability amplitude and cannot be interpreted in the 
same way as in the mean spectrum. 
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ity variations. Therefore, the fractional luminosity variations 
for SWIFT J1753. 5-0127 and GX 339-4 may be similar, but 
GX 339-4 shows a significantly greater normalisation in the 
PSDs shown in figure [T] 

4.2 Variable coronal geometry 

We have seen that it is likely that reprocessing of the power- 
law drives at least some of the disc variability seen in hard 
state sources, and possibly all of it on time-scales < 1 s. How- 
ever, model-independent covariance ratio plots and spectral 
fitting show that the covariance spectra of both sources 
demonstrate increased disc blackbody variability with re- 
spect to the power law on longer timescales. There are sev- 
eral possible explanations for this pattern, but the key thing 
that any successful model needs to achieve is an increase in 
disk variability on longer timescales without a concomitant 
rise in power law variability. Also, it is important to note 
that the additional disc variability must still be correlated 
with power-law variations, because the enhanced disc varia- 
tions appear in the covariance spectrum, which is identical to 
the rms spectrum over the high signal-to-noise energy range 
covered by the disc (i.e., coherence is unity). If the disc vari- 
ations were independent of the power-law they would cancel 
to some extent, since they would be uncorrelated with the 
power- law component, and covariance would be smaller than 
the rms. Thus the blackbody variations map on to power-law 
variations but with larger amplitude. 

One possibility is to change the geometry of the system 
on longer time-scales. For example, a variable coronal scale 
height on longer timescales could lead to changes in the solid 
angle of disc heated by the power-law, thus increasing the 
disc blackbody variability amplitude relative to the power- 
law. Weaker correlated power-law variation could then be 
produced if scale-height correlates with power-law luminos- 
ity. Alternatively, correlated power-law variations could be 
due to variable seed photon numbers from the disc due to 
the enhanced variable heating, but in either case the addi- 
tional power-law variability must be of smaller amplitude 
than the observed blackbody variability. Regardless of these 
model-dependent arguments, one can make a simple obser- 
vational test of the variable-geometry model, by comparing 
the variability of reflection on long and short time-scales. 

GX 339-4 shows significant reflection features in its 
covariance spectra, so any variation in coronal geometry 
should manifest itself as increased variability in the reflec- 
tion components as the disk sees the varying power law. 
However, the spectral fit parameters given in Table [2] indi- 
cate that there is little change in reflection amplitude be- 
tween long and short time-scales, both in terms of reflection 
covering fraction and iron line equivalent width (e.g. ratio of 
line flux to power-law normalisation, which is meaningful be- 
cause the power- law indices are so similar). To highlight this 
similarity in the covariance spectra, we show in figure [7] the 
data/model ratios for the PCA spectra with the reflection 
components taken out of the fit. The spectra demonstrate 
no significant change in the reflection continuum and asso- 
ciated iron line. This implies that if there is any change in 
geometry, it is small, and the increase in the thermal com- 
ponent of the variability spectra on longer timescales has a 
different origin. 

To place this result on a more rigorous statistical foot- 
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Figure 7. Data to model ratios of covariance spectra in GX 339-4 
with reflection normalisations set to zero 

ing, we show in figure [8] a contour plot of the best-fitting 
short and long-time-scale reflection covering fractions, which 
are obtained only from fits to the PCA data, which are most 
sensitive to the reflection continuum^ If variable geometry 
is the cause of the enhanced long-term blackbody variabil- 
ity, we would expect the covering fraction to show a similar 
enhancement in variability on long time-scales. The dashed 
line in the figure shows the largest ratio of short-to-long- 
time-scale covering fraction which crosses the 99 per cent 
confidence contour, with a value of 0.81, placing a 99 per 
cent confidence upper limit of 23 per cent on any increase 
in the covering fraction on long time-scales compared to 
short time-scales (the 90 per cent confidence upper limit 
on any increase is 12 per cent). In contrast, the long-time- 
scale blackbody variability amplitude increases by 30 per 
cent compared to that on shorter time-scales, which can- 
not be explained by the permitted increase in variable re- 
flection at greater than the 99 per cent confidence level. 
Thus, in GX 339-4 at least, we can rule out long-time-scale 
changes in coronal geometry as a viable explanation for the 
enhanced blackbody variability amplitude. Note also that 
the same arguments apply to any other non-geometric argu- 
ments which seek to produce the extra long-term disc vari- 
ability by varying the power-law contribution as seen by the 
disc, e.g through variable beaming of the power-law towards 
the disc. Due to the weakness of the iron line in SWIFT 
J1753. 5-0127 it is not possible to place similar constraints 
on coronal geometry changes in this source, although by 
analogy we expect the same interpretation to apply. 

4.3 Intrinsic disc variability 

Having established that variable coronal geometry cannot 
explain the extra blackbody variability on longer timescales, 
we are left to consider the possibility that the variations 

9 As a result, the best-fitting covering fraction is significantly 
lower than observed in the joint fit with EPIC-pn, but since we 
are interested in the relative change in reflection between different 
time-scales, this does not affect our basic result. 
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Figure 8. GX 339-4 confidence contour plots for the permit- 
ted range of reflection covering fraction in simultaneous fits to 
the long and short-time-scale covariance spectra. The dotted line 
passes through the 99 per cent confidence contour, so its gradient 
defines an upper limit to the ratio of long to short-time-scale cov- 
ering fraction. Also shown are the 90 and 68 per cent confidence 
level contours. 

are intrinsic to the disc itself. Perhaps the simplest pos- 
sibility is that the fluctuations are due to accretion rate 
fluctuations in the disc which undergo viscous damping be- 
fore they reach t he corona. Such damp ing is expected in 
thin discs (e.g. see lChurazov et aDl200l| ). but will not be as 
significant in geometrically thick flows, which might corre- 
spond to the corona. Therefore one can envisage that the 
long time-scale variability, e.g. corresponding to the low- 
frequency Lorentzian in the PSD, is generated in the thin 
disc, producing relatively large amplitude blackbody vari- 
ations but being damped before reaching an inner coronal 
emitting region (perhaps inside the innermost radius of the 
thin disc). The correlated power-law variability could then 
arise either from the residual undamped variations in ac- 
cretion rate which reach the corona, or be driven by seed 
photon variations from the disc. 

If the inner radius of the disc were to fluctuate on long 
time-scales, then the varying disc area would also introduce 
extra blackbody variability. This model is in some sense 
analogous to that of changes in coronal scale height, since 
any change in disc area will vary the solid angle of disc seen 
by the corona which will cause an increase in reflection as 
well as correlated power-law variability driven by seed pho- 
ton variations. One could mitigate these effects if the coro- 
nal properties were linked to those of the disc inner-radius, 
e.g. as inner radius decreases, so does coronal scale- height. 
This situation might be expected if the corona is formed b y 
evaporation of the disc (|Mever-Hofmeister fc Meveij|20o3) . 
so that condensation of the disc will drain and cool the coro- 
nal plasma, leading to a reduction in scale-height. 

4.4 Wider implications 

We have established that the hard-state variability on 
time-scales greater than seconds, corresponding to the low- 
frequency Lorentzian PSD component, is very probably pro- 
duced by variations intrinsic to the accretion disc, perhaps 
in the form of propagating accretion rate fluctuations, as 



envisag ed bv lLvubarskiil (|l997l ) to explain the broad shapes 
of observed PSDs. These variations then manifest as weaker 
power-law variations, either thro ugh propagation of ac cre- 
tion fluctuations to the corona |Ch urazov et al.l l200ll ) , or 
through variations of seed photons from the disc, which are 
Compton upscattered in the corona. On shorter time-scales, 
the blackbody variations are probably mostly produced by 
X-ray heating of the disc by the power-law, which is a re- 
quired outcome of the X-ray reflection directly observed in 
these systems. The exact split between direct emission from 
intrinsic disc variations and X-ray heating is difficult to 
judge, since some residual intrinsic variations may remain 
on short time-scales, and X-ray heating will also contribute 
on long time-scales. But in the case of intrinsic disc accretion 
fluctuations, the variable X-ray heating is itself a product of 
those fluctuations, and so any intrinsic disc variations on 
long time-scales must be large, at least comparable to the 
fractional rms of blackbody emission, i.e., 40 per cent in 
both sources. 

It is interesting to contrast the large intrinsic disc vari- 
ability in these hard state sources with that in the soft states, 
where the disc emission dominates the bolometric luminos- 
ity. The soft sta tes are well-known for showing very weak, if 
any, variability dHoman et al.|[200ll ). and the strongest vari- 
ability which is seen, e.g. in Cyg X-l, is associated with 
the power-law, with the disc being remarkably constant 
l|Churazov et al.ll200ll ). Therefore it seems likely that hard 
state discs are inherently unstable compared to soft state 
discs. This difference may represent just another observable 
distinction between hard and soft states, but it is interesting 
to speculate that it may play a more primary role in creat- 
ing the other observed differences, such as a strong corona 
and jet formation in the hard state. Certainly, it seems likely 
that intrinsic disc variability plays an important role in de- 
termining the PSD shape in the hard state. For instance, 
the low energy PS D of GRS 1915+105 in the x-class hard 
intermediate state ^Rodriguez et al] 120041 ) shows a low fre- 
quency component which disappears at higher energies, pos- 
sibly indicating that this source is demonstrating intrinsic 
disc variability. 

It is worth noting here that X-ray/optical studies of 
AGN also show evidence for reprocessing on short time- 
scales and intrinsic disc variations on longer time-scale s 
jUttlev et al.ll2003l ; lArivalo et al.ll2008l ; iBreedt et al.ll2009T ). 
However these AGN are relatively luminous and radio-quiet 
and so are likely analogues of BHXRB soft states, perhaps 
indicating that disc-stability shows a mass-dependence, e.g. 
related to the tra nsition between gas- and radiation- pressure 
dominated discs l)Maccarone. Gallo fc Fenderll2003h . 

Our results strongly suggest that disc variations are 
responsible for the low-frequency component in the hard 
state PSD. A number of authors have suggested that the 
low-frequency Lorentzian corresponds to the viscous time- 
scale of the inner, truncation radius of the thin disc (e.g. 
IChurazov et al.ll200ll : lDone et al.ll2007f >. a view which is con- 
sistent with our results. The higher-frequency PSD compo- 
nents may then be produced in the corona, which is likely 
to be geometrically thick so will show naturally shorter vari- 
ability time-scales. The viscous time-scale scales with scale- 
height (H) over radius (R) as (H/R) 2 . Assuming the ratio 
of disc scale-height to radius H/R ~ a ~ 0.1 (where a is the 
viscosity parameter), the predicted inner disc radius corre- 
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sponding to the observed low-frequency PS D peak around 
~ 0.05 Hz is 15 R G (see iDone et all 120071 '). However, this 
radius may be even smaller for smaller H/R, which would 
make the disc truncation radius consistent with the results 
from fits to the iron line (this work, and iReis et al.l [20081 . 
2009). However, in the latter case the corona would need 
to be very compact, or seed photon variations from intrin- 
sic disc variability would modulate the power-law with a 
similar amplitude to the disc. In either case, assuming the 
disc is thin and that it varies on the viscous time-scale, 
its inner radius must be relatively small. Such a picture 
is very different from earlier models for the hard state, 
where the disc is highly truncated and power-law emis- 
sion is produced by a very extende d corona or ADAF (e.g. 
lEsin. McClintock, fc NaravarJll997r i. 



disc on viscous timescales, is that the disc truncation ra- 
dius is < 20 Rg- 

This work highlights the importance of measuring spec- 
tral variability on a range of time-scales. Mean spectra, 
which describe the average properties of a source, provide 
no information on how different spectral components are re- 
lated to one another as a function of time. By using the 
covariance spectra we have been able to disentangle the cor- 
related spectral components in these two sources, identify 
thermal reprocessing as the mechanism by which variability 
is correlated in different bands, produce model-independent 
evidence for additional blackbody variability on longer time- 
scales and therefore associate intrinsic disc variability with 
the low frequency Lorentzian feature seen in hard-state 
PSDs. 



5 CONCLUSIONS 

Spectral fits to observations of BHXRB sources in the hard 
state show increasing evidence for both power law and black 
body components. In this work we have explored the hard- 
state variability of two sources with known soft excesses, 
SWIFT J1753.5-0127 and GX 339-4. Our findings are sum- 
marised below. 

1. We have introduced a new spectral analysis technique, 
the covariance spectrum, which measures the correlated 
variability in different energy bands. This technique over- 
comes the problems of low signal-to-noise and bias asso- 
ciated with the rms spectrum and has smaller statistical 
errors. 

2. PSDs of the two sources demonstrate larger low- 
frequency power in the soft band. 

3. The longer time-scale (2.7-270 s) covariance spectra of 
both sources are softer than the short time-scale (0.1-4 s) 
covariance spectra, due to additional disc variability, i.e., 
extra disc variability occurs on longer time-scales with- 
out a concomitant rise in power law variability on such 
timescales. However, the coherence of the rms and co- 
variance spectra show clearly that disc variations are not 
independent of the power law variations. 

4. The strength of reflection features that are detected 
in the short time-scale covariance spectra of GX 339-4 
are consistent with the observed blackbody variations on 
those time-scales being driven by thermal reprocessing of 
the power-law emission absorbed by the disc. However, 
the reflection covering fraction and iron line equivalent 
width show little change between short and long time- 
scales, implying that additional reprocessing, due to coro- 
nal geometry change, is not responsible for the additional 
blackbody variability seen on longer time-scales. 

5. The extra blackbody variability seen on longer time- 
scales appears to be intrinsic to the accretion disc it- 
self, giving rise to the extra low-frequency power in the 
PSD. This represents the first clear evidence that the 
low-frequency Lorentzian component in hard state PSDs 
is produced by disc variability. Models invoking damped 
mass accretion rate variations or oscillations in the disc 
truncation radius can satisfactorily explain the observed 
pattern of variability. 

6. The implication of such variations occurring in a thin 
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